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ABSTRACT 



So far, solar-like oscillations have been studied using ra- 
dial velocity and/or light curve variations, which reveal 
frequencies of the oscillation modes. Line-profile vari- 
ations, however, are also a valuable diagnostic to char- 
acterise radial and non-radial oscillations, including fre- 
quencies, amplitudes, the spherical mode wavenumbers 
(£, to) and the stellar inclination angle. Here we present 
a line profile analysis of e Ophiuchi, which is a pulsat- 
ing red giant. The main differences compared to previous 
line profile analyses done for heat-driven oscillations are 
the small amplitudes and the predicted short damping and 
re-excitation times in red giants. 

Two line diagnostics have been tested to see whether 
these are sensitive to the small line profile variations 
present in red giants. In addition, line profiles have been 
simulated with short damping and re-excitation times and 
are compared with the observations. This comparison re- 
veals that non-radial modes are detected in the observed 
line profile variations of e Ophiuchi. This is rather sur- 
prising, as theoretical predictions favours the occurrence 
of radial modes. 
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1. INTRODUCTION 



The discovery of oscillations in the sun, about five 
decades ago, was the onset of seismology of stars. Os- 
cillations probe the interior of the stars and are therefore 
a direct means to determine their internal structure. Solar- 
like oscillations are excited by turbulent convection near 
the surface of cool stars of spectral type F, G, K or M. 
They show radial velocity variations with amplitudes of 
typically a few cms -1 to a few ms _1 , and with periods 
ranging from a few minutes for main sequence stars to 
about half an hour for subgiants and a couple of hours for 
giants. 



Table 1. Basic stellar parameters of c Ophiuchi taken 
from Jll/: Effective temperature (T e g) in Kelvin, rota- 
tional velocity (vshiij in kms~ x , parallax (ir) in mas, 
distance (d) in pc, the apparent magnitude (m v ) and ab- 
solute magnitude (My) in the V 'band. 



parameter 


e Ophiuchi 


T eff [K] 


4887 ± 100 


usini [kms -1 ] 


3.4 ±0.5 


7r [mas] 


30.34 ±0.79 


d [pc] 


33.0 ±0.9 


m v [mag] 


3.24 ±0.02 


My [mag] 


0.65 ±0.06 



During the last decade, the refinement of the techniques 
to perform very accurate radial velocity observations, 
made it also possible to observe solar-like oscillations in 
distant stars. Recently, these solar-like oscillations are 
indeed observed in several red giants, for example in £ 
Hydrae JU, e Ophiuchi |2] and r\ Serpentis @|. 

This proceedingspaper describes observations as well as 
a frequency analysis of different parameters and a line 
profile analysis of e Ophiuchi. It is organised as follows. 
In Section 2, we start with a discussion on the target se- 
lection, followed in Section 3 by the description of the 
observations. In Section 4 the data analysis is described, 
i.e. first, the frequency determination of different oscilla- 
tion diagnostics and subsequently the mode identification 
using these frequencies. In Section 5 we show simula- 
tions, including the short damping and re-excitation times 
present in solar-like oscillations of red giants |4]. Finally, 
in Section 6 we draw some conclusions. 



2. TARGET SELECTION 



In order to obtain radial velocity variations with an ac- 
curacy of order ms _1 , the target selection as well as the 
observing strategy are important. First, the stars have to 
be bright enough to obtain high signal to noise ratios to 
detect the signal. Secondly, the rotational velocities have 



to be low in order to have narrow spectral lines, which 
are needed for precise radial velocity measurement. Also 
other infering phenomena like star spots or companions 
should be avoided. 

For the observation strategy we have to take into account 
that the expected (and observed) period of the solar-like 
oscillations in red giants is of order a couple of hours. 
This means that the observation times should be of order 
minutes, in order not to average over a large part of the 
oscillation cycle. Furthermore single site observations al- 
ways introduce one day aliases due to the diurnal cycle. 
Multi site observations, or continuous observations from 
space can significantly reduce this aliasing effect. 

The G9.5III giant e Ophiuchi was selected for observa- 
tions of solar-like oscillations, because of its brightness, 
m v = 3.24 mag, its low rotational velocity, vsvai = 3.4 
km s~\ and the low declination, 5 = -04 41 33.0 |5]. 
The latter was very important to be able to observe the 
star from both La Silla, Chile and Observatoire de Haute 
Provence, France. More basic stellar parameters of e 
Ophiuchi are shown in Tabled Furthermore the Hippar- 
cos catalog |6] gives a photometric variation of less than 
0.06 mag and the star is mentioned in the '1988 revised 
MK spectral standards for stars GO and later' |7]. There- 
fore the presence of starspots or companions is not very 
likely. 



3. OBSERVATIONS 



The observations of e Ophiuchi were performed in a 
bi-site campaign with the fibre-fed echelle spectrograph 
CORALIE, mounted on the Swiss 1.2 m Euler telescope 
at La Silla (ESO, Chile), and the fibre-fed echelle spectro- 
graph ELODIE, mounted on the French 1 .93 m telescope 
at the Observatoire de Haute Provence (France), during 
the summer of 2003. The spectra range from 387.5 nm to 
682 nm in wavelength. The observation times were kept 
short enough, in order not to average over a large part of 
the oscillation, but long enough to reach a signal to noise 
ratio of at least 100 at 550 nm. 

The radial velocities of e Ophiuchi were determined with 
the optimum weight method 1 8 ] and are available in 0] • 
Variations of order 20 m s _1 due to stellar oscillations 
were clearly detected by the authors. We refer the reader 
to 1 2] for more information on the radial velocity analysis 
of the star. 

In this work, we use the same data as in 0], but we fo- 
cus on the line profile variations of the star rather than on 
the radial velocity. Such variations can be derived with 
high accuracy from a single unblended line provided that 
it is not subjected to Stark broadening (see, e.g., |9] and 
references therein) and that the spectra have a high signal- 
to-noise ratio. The latter condition is not fulfilled for our 
data of e Ophiuchi. Therefore, we perform a line profile 
analysis from a cross correlation profile computed of each 
spectrum. This was done by using a box-shaped mask 
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Figure 1. Power spectrum of e Ophiuchi. The lower black 
one is obtained by [2] from radial velocities derived with 
the optimum weight method [S], the grey one in the mid- 
dle is obtained from (v) and the top one is obtained from 
the bisector velocity span. For clarity the latter two are 
shifted. The dotted vertical lines indicate the dominant 
frequencies in (v). These are used in Figure^\to compute 
the amplitude across the profile. 



11011 . This cross correlation profile has an appreciably 
higher signal-to-noise ratio compared to a single spectral 
line. Moreover, Willi have shown that a line profile analy- 
sis in terms of stellar oscillations is possible using a cross 
correlation profile rather than individual lines. This is a 
similar conclusion to the one by 1131 . who showed that 
the cross-correlation bisector contains the same informa- 
tion as single-line bisectors. 



4. DATA ANALYSIS 



As mentioned earlier, |2|] used the optimum weight 
method |8] to calculate the radial velocities. This is in- 
trinsically the most precise method, because the whole 
spectrum is used. With this method the observation with 
the highest signal to noise ratio of each night is used as a 
reference. These nightly reference points are effectively 
a high pass filter and therefore the low frequencies in the 
power spectrum drop to a level near zero. The resulting 
power spectrum of e Ophiuchi is shown as the lower black 
graph in Figure^ 

From this power spectrum, [2] determined two possible 
large separations: 4.8 /iHz or 6.7 /iHz. Two different stel- 
lar evolution models fit these large separations. A stellar 
model with a mixing length of convection a = 1.6 (ex- 
pressed in local pressure scale heights) and a mass of ap- 
proximately 2.8 solar masses has a large separation close 
to 6.7 /iHz, assuming that the detected modes are radial. 
A second stellar model with a higher value of a = 1.8 
and a lower mass of approximately 1 .9 solar masses has 
a large separation of 4.8 /iHz. The data in [2] did not al- 
low to discriminate between these two values of the large 
separations from theoretical arguments alone. 



Recently, e Ophiuchi was also observed with the MOST 
satellite. These continuous photometric observations, de- 
scribed elsewhere in this volume [3 , rule out the 6.7 
fiHz large separation and thus reveal the evolutionary 
phase of e Ophiuchi, again assuming that its detected 
modes are all radial. 



4.1. Moments 



Although the radial velocities derived from the optimum 
weight method are intrinsically the most accurate ones, 
this method is not useful for mode identification. In or- 
der to get the best line profile analysis a cross correlation 
profile of each spectrum is calculated, using a box-shaped 
mask 1 10]. This cross correlation profile has an increased 
signal to noise ratio compared to a single spectral line. 
Such a profile can be described by its moments and one 
usually considers the first three moments for mode iden- 
tification 1 14]. The first moment (v) represents the cen- 
troid velocity of the line profile, the second moment (v 2 ) 
the width of the line profile and the third moment (v 3 ) the 
skewness of the line profile. The quantity (v) is a partic- 
ular measure of the radial velocity and should therefore 
show similar frequency behaviour as the radial velocity 
derived from the optimum weight method. The power 
spectrum of (v), where the average value of (v) per night 
is used as a reference, is shown in the middle of Figure^ 
Despite the differences in the methods to compute a mea- 
sure of the radial velocity, both power spectra reveal the 
same dominant frequencies. 



4.2. Bisector Velocity Span 



4.3. Mode identification 

In order to characterise the wavenumbers (I, m) of the 
oscillation modes, the amplitudes across the line profile 
are compared with the ones obtained from simulations. 
These amplitudes are determined by fitting a harmonic 
function, with the dominant frequency of (v), to the flux 
values at each velocity pixel of the time series of spec- 
tra fl6ll . This is shown schematically in Figure |2] for 
{I, m) = (0, 0), and in FigureEJfor (£, m) = (2, 2). The 
left panel shows profiles obtained at different times with 
an arbitrary flux shift. The dotted / dashed lines indi- 
cate two examples of velocity points at which a harmonic 
function is fitted. A harmonic fit at the centre of the pro- 
files is shown in the top middle panel of Figures Eland|3] as 
a function of phase. A harmonic fit through a wing of the 
profiles is shown in the bottom middle panel of Figures^ 
and [3] A difference between the amplitudes of the har- 
monic fit at the different velocity points is clearly seen. 
The behaviour of the amplitude across the whole profile 
is shown in the top right panel, while the phase across the 
whole profile is shown in the bottom right panel. This 
diagnostic is suitable to be compared with a similar one 
derived from the data. 

The observed amplitude and phase distributions of e 
Ophiuchi for its three dominant frequencies are shown in 
Figure|4] as a function of velocity across the line profile. 
If only radial modes would be dominant in e Ophiuchi's 
line profile variations, then all three amplitude distribu- 
tions should have the same shape. However, the detected 
shapes of these amplitude diagrams clearly differ for the 
three frequencies. Moreover, the very characteristic am- 
plitude shape for an axisymmetric (m = 0) mode is not 
recovered for the dominant mode of the star. This im- 
plies that at least one non-radial mode is detected in the 
line profiles of e Ophiuchi. 



A line bisector is a measure of the displacement of the 
centre of the red and blue wing from the core of the spec- 
tral line at each residual flux. The bisector velocity span 
is defined as the horizontal distance between the bisector 
positions at fractional flux levels in the top and bottom 
part of the line profile, see for instance 1 1 5 ] . 

The bisector velocity span is a measure of the line profile, 
mostly used as a tool to distinguish between companions 
and intrinsic stellar activity as sources of the observed 
radial velocity variations. The Doppler shift induced by a 
companion shifts the entire spectrum, but does not change 
the spectral line shapes, while oscillations and star spots 
do change the spectral line shapes. 

No dominant frequencies are obtained from the power 
spectrum of the bisector velocity span, which is plotted 
as the top graph in Figure [2 We conclude that for oscil- 
lations with low amplitudes, the bisector velocity span is 
not a useful diagnostic to distinguish between planetary 
companions and intrinsic variations in a star. This result 
is consistent with the results recently found by 1 12]. 



5. SIMULATIONS OF LINE PROFILES 



The results obtained in the previous section are not an- 
ticipated from a theoretical viewpoint (e.g. 1 17]). One 
may therefore wonder if the profiles simulated for infi- 
nite mode lifetimes differ too much from reality, i.e. from 
those for finite lifetimes as in the data. The damping and 
re-excitation times for red giant oscillations turn out to 
be short. Indeed, |4] derived an oscillation mode lifetime 
for £ Hydrae (G7III) of only approximately two days. 
In order to take this into account, and to investigate the 
robustness of our conclusion against finite lifetimes, we 
computed line profile variations simulated for stochasti- 
cally excited modes. Such line profile variations are not 
yet available in the literature. We provide them here in an 
attempt to make a definite conclusion on the nature of e 
Ophiuchi's dominant oscillation modes. 

A damped and re-excited oscillation mode is damped by 
a factor e -??t , with rj the damping rate, and re-excited 
before it is able to damp out. As a consequence both 




Figure 2. Schematic representation of the amplitude distribution across the profile for simulated data with (£,m) — (0,0), 
an amplitude of the pulsation velocity of 0.04 Ions -1 , an inclination angle of 35°, and intrinsic line width of 4 kms^ 1 
and a vsini of 3.5 kms . Left: Profiles obtained at different times are shown with an arbitrary flux shift. The dashed 
and dotted lines indicate the two velocity values at which the harmonic fits shown in the two middle panels are obtained. 
Middle top: harmonic fit at the centre of the profiles. Middle bottom: harmonic fit at a wing of the profiles. Right top: 
amplitude across the whole profile. Right bottom: phase across the whole profile. 




Figure 3. Schematic representation of the amplitude distribution across the profile for simulated data with (£, m) — (2, 2), 
an amplitude of the pulsation velocity of 0.04 kms -1 , an inclination angle of 35°, and intrinsic line width of 4 kms^ 1 
and a vsini of 3.5 kms . Left: Profiles obtained at different times are shown with an arbitrary flux shift. The dashed 
and dotted lines indicate the two velocity values at which the harmonic fits shown in the two middle panels are obtained. 
Middle top: harmonic fit at the centre of the profiles. Middle bottom: harmonic fit at a wing of the profiles. Right top: 
amplitude across the whole profile. Right bottom: phase across the whole profile. 
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Figure 4. Amplitude (top) and phase (bottom) distributions (thick lines) and their uncertainty (indicated in grey) as a 
function of velocity across the line profile of e Ophiuchi for three significant frequencies obtained from (v); u/ v \ — 5.03 
c/d(58.2/j,Hz) (left) V( v ) =5.46 c/d(63.2fiHz) (middle) and =5.83 c/d (67. &fiHz) (right). The mean radial velocity 
of the star, is found to be approximately —9.4 km/s. 
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Figure 5. Amplitude and phase distributions for simulated line profiles with a two day damping time at an inclination 
angle of i — 55°. For each mode, ten different realisations are shown in grey and the average distribution is shown in 
black. Top left: i = 0, m = 0, top right: I = 1, m = 0, middle left: i = 1, in = 1, middle right: £ = 2, m = 0, bottom 
left: I = 2, m = 1, bottom right: I — 2, m = 2. 




Figure 6. Amplitude and phase distributions for simulated line profiles with a two day damping time at an inclination 
angle of i = 75°. For each mode, ten different realisations are shown in grey and the average distribution is shown in 
black. Top left: £ — 0, m — 0, top right: I = 1, m = 0, middle left: I = 1, m = 1, middle right: I — 2, m — 0, bottom 
left: I = 2, to = 1, bottom right: I = 2, m = 2. 



the amplitude and the phase of the oscillation are time 
dependent. 

/0) = A(t) sin(27rz4 + ^>(i)). (1) 

To simulate such an oscillator we follow the description 
of 1 18], and we compute 

f(t) = B(t) sin(27ri/£) + C(t) cos(2tt^), (2) 

where we let the amplitudes B and C vary with a first or- 
der autoregressive process in a discrete time domain with 
time step At: 

£ n = e-" At B n _i +e„+i, (3) 

where e n +i is a Gaussian distributed excitation kick. For 
more details we refer to I18ll . The above is applied to the 
equations for the pulsation velocity which are obtained 
by taking the time derivative of the displacement compo- 
nents mentioned in Section 3.2, Eq.(2) of 1 19]. In case 
rotation is neglected, the three spherical components of 
the pulsation velocity, i.e. v r , vq and v v , for the damped 
and re-excited case become: 

Vr = -V P NPPP (COS 8) e -n(t-nAt ktck ) # 

{B n s'm(mip + 2-Kvt) + C n cos(nrnp + 2m>t)), (4) 

v e = -^ Up ^ m ^(P™(cos0))e-"(*-" Atfc " fc ) * 

(B n s'm(mip + 2nvt) + C n cos(mip + 2irvt)) 1 (5) 

v v = -mKvpNF^PP (cos 6») e -"(*- nAtMc ^ * 
v y sm9 

(B n cos(nup + 2irvt) — C n sin(m(p + 2m/t)), (6) 



with v p proportional to the pulsation amplitude, N™ 
the normalisation factor for the spherical harmonics 
Y e m (6, f) = Pp (cos 6)e lm{ e and K the ratio of the hori- 
zontal to the vertical velocity amplitude. 

Line profiles of oscillations with different wavenumbers 
(£, to), including damping and re-excitation, are simu- 
lated at the observation times of e Ophiuchi and according 
to its observed amplitudes. Subsequently, (v) is deter- 
mined for these profiles in the same way as for the obser- 
vations. The amplitude and phase distributions are then 
computed for the dominant oscillation frequencies, as de- 
scribed in section 4.3. 

Only in case of an infinite number of observations over 
an infinite timespan, the real dominant frequency can be 
obtained for a damped and re-excited oscillator. Smaller 
samples of observations, single realisations, are all dif- 
ferent and have their own dominant frequency and am- 
plitude distribution. In Figures |5] and |6] simulated ampli- 
tude and phase distributions for different wavenumbers 
are plotted at inclination angles i = 55° and i = 75°, 
respectively. For each mode, ten single realisations are 
plotted in grey, with the average plotted in black. From 
these Figures it becomes clear that, although all realisa- 
tions are different from each other, it is still possible from 
the amplitude distributions, to distinguish between differ- 
ent to values, just by visual inspection. In particular, the 
different realisations for to = modes all lead to zero 
amplitude in the line centre. In contrast with the pro- 
nounced shapes of the phase diagrams of oscillations with 
infinite lifetimes (see e.g. right bottom panel Figures |2] 
and|3]or Il20l0 . Figures |5] and |6] show that the phase dis- 



tributions across a profile can be quite different for each 
realisation due to damping and re-excitation. Therefore, 
phase distributions across a profile of oscillations with fi- 
nite mode lifetimes can not yet be explored in terms of 
mode identification. These results are robust against a 
change in inclination angle, as can be seen by comparing 
Figures |5]and[6] with a similar plot for a different inclina- 
tion provided in ll20ll . 



6. CONCLUSIONS 



We have found a different amplitude distribution across 
the cross correlation profile for the dominant frequencies 
of the pulsating red giant star e Ophiuchi. This clearly 
indicates that modes with different wavenumbers (£, m) 
must be present in the data. Moreover, the amplitude dis- 
tribution of the dominant mode is clearly different from 
the one due to m — modes. Thus, we have detected 
non-radial modes in the observed cross correlation pro- 
files of this star. 

We compared the observed amplitude distributions of e 
Ophiuchi with the one of simulated cross correlation pro- 
files in an attempt to identify the wavenumbers of the in- 
dividual modes. Due to the fact that single realisations 
of different modes are in some cases comparable, e.g. for 
(£, m) = (1,1) and (2,1) (see Figures[5] and [6ji, differ- 
ent sets of wavenumbers are equally likely for observed 
modes. It is therefore clearly not yet possible to iden- 
tify the individual modes unambiguously with the method 
used in the present work, but clear evidence for modes 
with ra/0 values is found. 

12011 also analysed the line profile variations of three ad- 
ditional pulsating giants and came to similar conclusion 
than for e Ophiuchi for two of these three stars. More- 
over, non-radial mixed oscillation modes were recently 
also found in a red giant observed with the MOST space 
mission (|21J, these proceedings). It therefore seems 
that red giant stars pulsate in both radial and non-radial 
modes. These observational results constitute a fruit- 
ful starting point to fine-tune theoretical predictions of 
stochastically excited modes in evolved stars. 
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